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^ ' Abstract 

a^ 

_^ ' It has been shown that the high-temperature hmit of perturbative thermal QCD 

is easily obtained from the Boltzmann transport equation for 'classical' coloured 
Qh! particles [2]. We generalize this treatment to curved space-time. We are thus able to 

construct the effective stress-energy tensor. We give a construction for an effective 
action. As an example of the convenience of the Boltzmann method, we derive the 
high-temperature 3-graviton function. We discuss the static case. 



1 Introduction 



The high-temperature behaviour of one-loop diagrams in perturbative thermal 
QCD is T^ (where T is the temperature), for any number of external gluon 
lines (we ignore diagrams with external quark lines in this paper) [1]. An 
elegant recent paper [2] has shown how the colour current Jl^{x;A) (where 
y4"(x) is the gluon field) which corresponds to this infinite set of graphs may 
be obtained from the Boltzmann equation for classical coloured particles [3] 
moving in an external classical gluon field. It is also known [1,4,5] that there 
is a nonlocal action iy(y4) such that 

, SW 

y/^j^{x; A, g) = - ^ ^ . . . (1.1) 
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(We have inserted a metric g here anticipating step (i) below.) So far as we 
know, it is not obvious from the Boltzmann equation approach [2] that (1.1) 
should be integrable. 

Here we generalize the work of [2] in the following ways: 

(i) We generalize to curved spacetime, with metric gfj_^. This allows us to treat 
graphs with just external graviton lines (with T'^ behaviour), and also graphs 
with both gluon and graviton lines (T^ behaviour) 

(ii) As well as the current (1.1), we are then able to define an energy tensor 
by 

, 6W 

,/^T^-^(x;A,g) = -2—^^. (1.2) 



(iii) We write down an exact, but implicit, solution to the Boltzmann equation. 

(iv) We explain why equations (1.1) and (1.2) are integrable, and obtain an 
expression for W in terms of the stationary value of the (local) classical action 
for an ensemble of classical particle trajectories in the gluon and gravitational 
fields. This action contains both the T'^ and the T^ terms. 



(v) We discuss the positivity of T^*^. 

(vi) As an example of the effectiveness of the method, we derive the high- 
temperature 3-graviton function, which is quite complicated to obtain directly 
from perturbative quantum field theory. 

(vii) We discuss the static case, in which explicit solutions of the Boltzmann 
equation are easy to write down. 



2 The Boltzmann equation and the effective action 



Following [2], we imagine an ensemble of 'classical', coloured particles, follow- 
ing null trajectories x^{d), where 9 is an afiine parameter (which we take to 
have dimensions of length). Each particle has a 'classical' colour charge I°'{9) 
[3]. Let p^ = l~^x^, where / is a parameter with dimensions of length, which 
we use merely to keep track of dimensions. The distribution function f{x,p, I) 
in a generalized phase-space obeys the Boltzmann equation 



^"h - r»^|x -^"^;^(-)^"(-)^''|x 



+ r'"'iy^A\{x)^^ 



/ = o. 



(2.1) 



Here we have absorbed the QCD couphng constant into A (in order to 
avoid confusion with g = det((?^,y)). We also take the dimensions of A to 
be (length)"^ and 



^;. = d,K - d^Al + r'^A^Al. 



(2.2) 



Ahhough (2.1) is classical in the sense that it does not contain /i, it is necessary 
to remember that p does not have the dimensions of momentum: fi]) does. 

Equation (2.1) may be written in the "divergence" form 



^ ^--— r»^ 



d 



-p' 



dx^^ dp 



A t^^i 



dp- 
d 



rF^g^Y 



,A tJ-^^ 



_l_ fahcja 11 Ah 



g{x)f = 0. 



This is a sort of Liouville theorem, and shows that 
d'^xd'^pdlg{x) 



(2.3) 



(2.4) 



is an invariant phase-space element along a trajectory (for the definition of dl 
see [2]). The current and energy tensor are then given by 



fix) 



jd'pdi^/^p^rf{x,pj), 



T^^{x) = / d^pdl^/^p''p''f{x,p,l). 



(2.5) 



(2.6) 



From (2.1) follow the conservation equations 

D,r + r'^Air = o, 

D,T^^ + Frguxf = 0, 



(2.7) 
(2.8) 



where D denotes the coordinate (not colour) covariant derivative. 



Consider now a family of trajectories {x{9),p{9),I{9)), and define (Y,Q,Iq) 
by the conditions 

xr^Y, p^ r^Q, I ^ lo (as a; ^ oo) (2.9) 



and 



^^ = ^"^- (2-10) 



Thus three of the components of Y parametrize the position of the trajectory 
at infinity, and the fourth replaces 6 as an afiine parameter along a trajectory, 
y is a vector in the asymptotic Minkowski space (we assume that the colour 
and gravitational fields tend to zero at infinity). Then there are relations 

x = x{Y,Q,Io), p = p{Y,Q,Io), I = I{Y,Q,Io)- (2.11) 



Any function of a constant of the motion satisfies the Boltzmann equation 
(2.1), and an example of such a constant is Q" itself. We seek a solution 
of (2.1) which reduces to the equilibrium distribution at temperature T at 
infinity. The required solution is evidently 

f{x,p,I) = -^N[Q{x,p,I)] 2 P Sm 9{Qo), (2.12) 



where 



'^W = eMt^Q„/iT)±l P-1^' 



and C gives the number of spin and internal degrees of freedom. Here the 
factor h is because of the d'^p in (2.4), and the factor l^ is because of the l~^ 
in (2.9). Equation (2.12) is the exact solution of (2.1) with which we will work. 

We may now write equations for j^ and T'^'^: 



(fx^/^f^'{x)5A''^{x) = d^x / dV dl g{x)f{x,p, I) p^" P SA^^ix) 

d'Y / d^grf/oiV(Q)5(Q2)^(QO) 



(27r)3/2^ . 
X p'^iY, Q, /o)r(F, Q, Io)SA';ixiY, Q, Jo)). (2.14) 



So that 



^^nx) 



2Ch 
/2(27r)3 



where 

A(y,g,/o) 



det 



d''QdIoN{Q)5iQ')eiQ")A{Yix, Q, Jq), Q, Jo) 



X p^(F(x, g, Jo), g, /o)/o"(v'(x, g, /o), g, /o), (2.15) 



^ gx^(y,g,/o 



(2.16) 



Similarly 



-gT^'ix) 



2Ch 
/2(2^ 



d4grf/oiv(g),5(g2)^(gO)A(F(x, g, Jo), g, Jo 



xp>'^(F(x,g,/o),g,/o). 



(2.17) 



In order to evaluate these expressions in perturbation theory, one needs to 
find the functions in (2.12) in perturbation theory, and deduce perturbation 
expansions for Y{x, Q, Iq) and A. In contrast to the method of expansion 
used in [2], no derivatives of N{Q) appear in our method. (Of course, the two 
expansions must be equivalent when one uses the g-integration.) 

We now construct an effective action W{A,g) such that (2.15) and (2.17) 
are given by (1.1) and (1.2). There is no action for Wong's equations [3] in 
terms of /" (see for example [6] and references therein). Instead we introduce 
a quantity u, transforming say as the (complex) n-dimensional fundamental 
representation of SU(n) (which we now assume the colour group to be), with 
representation matrices T". This is related to I by 



m*T"m. 



(2.18) 



We consider -u to be a dynamical variable along the trajectory, with asymptotic 
value Mo- Then we take the action 



h 



d^Y 



S{Q,Uo;x,u;A,g) = 
{x^QfM + ilu*Qu} g^v{.x)x^x'' — ilu*u— 



lu*T''ux''A''^{x) 



(2.19) 



where the 'dot' differentiation was defined in (2.10). Thus x(Y) and u(Y) 
appear in (2.19) like 'fields' in y-space. The term in the brace in (2.19), 
being a differential, does not affect the equations of motion. The reason for 
introducing it will be explained later. In (2.19), factors of / are introduced to 



make everything dimensionless, and the factor h to give the whole thing the 
dimensions of an action. 

The equations of motion got by varying x{Y) and u{Y) in (2.19) are 

m + i:M;^(x)T"M = 0, (2.20) 

-d^-^)^^- ^^"^m(^) = 0' (2-21) 
where we have used (2.18). From (2.19) we deduce the Wong equation [3] 

and using this Eq. (2.21) gives the expected equation of motion. 
The next step is to define the effective action 

S{Q, Mo; A g) = S[Q, Mo; C(Q, ^o; A g),v{Q, Uo; A, g);A, g], (2.23) 

where x = ^ and u = v are solutions of the equations of motion (2.20) and 
(2.21) subject to the boundary conditions (2.9) and u — ^ Uq. Now make in- 
finitesimal variations in A and g. We have 

S{A + 5A,g + 5g) = S[^ + 6^,v + 6v; A + 5A,g + 6g] 

= S[^,v;A + 6A,g + 6g], (2.24) 

where 

^ + 6^ = ^{A + 6 A, g + 6g), v + 6v = v{A + 6 A, g + 6g), (2.25) 

and in the second equality in (2.24) we have made use of the fact that x = 
^, M = i; is a stationary value of S. Thus we deduce from (2.19) and (2.24) 
that 



^ I d'Y5Al{x{Y))rp^ = SiA + 6 A, g) - SiA, g) 

= S{A + 5A,g)-S{A,g). (2.26) 

Comparing this with (2.14), we see that (1.1) is satisfied by 

2C r 
W = -—- J d^Qduo6{Qy{Qo)N{Q)S{Q, uo; A, g). (2.27) 



An exactly similar argument shows that (2.27) also satisfies (1.2). 

We now explain the reason for inserting the brace into (2.19). The proof of the 
second equality in (2.24) involves an integration by parts, in which we require 
that 

ld'YQ. — [{g,,e - Q,}5^ + l^{v* - v*,}5v] = 0. (2.28) 

This would not be true without the contributions from the brace in (2.19), 
since 6C, and Sv do not tend to zero fast enough for large Y. In fact the whole 
contribution to (2.27) comes from the brace in (2.19). The other terms in 
(2.19) are zero for a solution of (2.20) and for a null-trajectory solution of 
(2.21). 

In order to use (2.27) to get the effective actions of hard-thermal loops, we 
need to expand (2.23) in powers of /. The leading term is of order /~*^. To 
this order, W is independent of the gluon field A. The 0(1) term in (2.21) 
has to be solved exactly to find this contribution to W. The effective action is 
that for hot matter in curved spacetime. It is, from (2.13), of order h{T/h)'^. 
The next nonvanishing terms in (2.23) are order l^^. To find them, we may 
solve (2.21) iteratively in powers of /, but (2.22) must be solved to all orders. 
The contribution from (2.27) is order h{T/h)'^. This is the effective action for 
gluons in curved spacetime. One might continue the expansion in powers of /, 
but for massless bosons the lower order terms may be infrared divergent. 

Finally, note that T*"^ from (2.17) is positive. It follows that the leading, /~^ 
term is positive. But we cannot say anything about the positivity of the /~^ 
term on its own. 



3 The 3-graviton hard thermal loop 



In practice, iterative solution of (2.1) seems to be an efficient way of calculating 
hard thermal contributions. As an example, we derive the 3-graviton function 
(keeping just the graviton terms in (2.1)). We have found it convenient to use 
the variable p^ = g^u{.x)p^, instead of p^. We define 

F{x^,p,,P) = f{x>^,g,,{x)p\P) (3.1) 



Then the factor g{x) is omitted from (2.4), and the relevant part of (2.1) 
becomes 

]9" (d^ + r>7^) F (^,P, I) = 0. (3.2) 



Using the well-known expression of the Christofell symbol F^^, and defining 
the graviton field 



/,"/3 = ^"/3 _ ^"/3^ 



(3.3) 



the Boltzmann equation can be written in the form 
p ■ d F {x,p,I) = LF {x, p, I) 



(3.4) 



where p ■ d = rj'^'^p^Pv and L is the differential operator 

d 



2 



d-y(l)°''^ } PaPiS 



dpj 



r'^iPadp+Ppda) 



(3.5) 



which is linear in the graviton field 0"^^. 

It is now easy to find recursively the solution of the Boltzmann equation, 
which is given by 



F(°) = 2Ch9 (po) 5 {ri^'p^pu) N (po) 



_p(") 



.P 



-a^J'""- 



(3J) 



where F*^"^ is of order 



The n-point graviton functions can be obtained by functional differentiation 
of the generating action using the relation 



6W 



Sg 



fj,i/ 



-9T, 



flU 



d'^p 

(2^ 



■p^PuF {x, p) . 



(3.7) 



From Eqs. (3.5) and (3.6) we can see that all F*^") have degree (—2) in 
p. Making the rescaling p -^ p/\v\ — Q) the n-graviton functions, in 
the momenta space, can be written in terms of dimensionless functions 



n^i.i...M„.„ (k\ ■■■,k^)=pl ^An^^^^6 [k' + ... + k^ 



^ //ii!^i---/i„i^„ y'^ 1 ■ ■ ■ ) "^ 1 vj 



Q2=o' 



(3.^ 



where the integration is over the directions of the 3-vector Q and 



cn 



IpI d|p| 



27r2 7 exp(|;ip|/T) 



C- 



n^T^ 



3or 



(3.9) 



is the energy density of the thermal particles. 

The functions ffMiui---fMnVn{k'^y '^ k"',Q) denote the coefficient 
of 0'^2i'2 ^f^2^ . . .(ptJ-ni^n ^f^Ti^ jj^ ^]-^g Integrand of (3.7), and can be obtained in 

a systematic way solving Eq. (3.6) for n = 1, 2, ■ ■ ■. Let us first consider the 
two-graviton function. In momentum space, Eq. (3.6) gives 






0°^ (A;) 



\ P ■ k I p ■ k 



(3.10) 



where k ■ dp = kad/dpa- Using the on-shell constraint, we can see that the 
surface term in the last equation can be neglected because in the limit po ~^ oo 
the function A^ (po) in the first equation of (3.6) goes much faster to zero than 
Po ~ IpP- This will also be true for all n-graviton functions. From Eq. (3.10) 
one can readily obtain f^uai3 {k, Q)- Then, Eq. (3.8) yields the following result 
for the 2-graviton function 



n 



fiuafS 



(k) 



I^U'-^^ Q.k 



(3.11) 



Q2=o 



This result is in agreement with the hard thermal loop 2-graviton function 
obtained by standard Feynman diagrammatic calculation [7] . The term inside 
the round bracket is identical to the forward scattering amplitude of a hard 
on-shell thermal particle off two external graviton fields. 

Let us now consider the 3-graviton function. Using Eq. (3.6) with n = 2 and 
adding total derivatives we obtain in the momentum space 



p^p.i^^'^ = r^ (A;') r" {k') Uuapp. (k\ e p 

0a/3 ^^2) ^pa (^3) 



F(°)P • a 



+ 



VpVuPpPa [Pakl + Pfikl 

p ■ k^p ■ {k^ + k'^) 



p ■ k^ 



r-a 



PpPuPaP/3 

p-{k^ + A;3) 



(3.12) 



Making p —^ Q in the coefficient of 0"^^ (A;^) 0'"^ (k^) and using Eq. (3.8), we 
obtain 



where 



A,.^p,. [k,k,Q)--k-dQ\ ——k ■ Oq 



+ 



Q-k^ "^g-(A;2 + A;3) 

,3 

■/3 



QfiQuQpQcr [Qaka + Qf)k^ 



Q2=o 



g ■ A;3 g . (P + A;3) 
+ (A;2, a, ^) ^-^ (P, p, a) . (3.14) 



A lengthy calculation (we have used computer algebra) has shown explicitly 
that (3.14) is the same as what is obtained from high-temperature perturbative 
quantum field theory and is therefore symmetric under all permutations of 

{k\^i,v), {k',a,(3), {k\p,a). (3.15) 



This verification makes use of eikonal identities like 

111 , , 

+ ^ ..^ ... =0- (3-16) 



Q-k^Q-k^ Q-k^Q-k^ Q-k^Q-k^ 

We stress that it would be very hard to guess, directly from perturbative 
quantum field theory that the 3-graviton function could be written in such a 
simple form as given by (3.14) (compare with the second work in reference 

[7]). 



4 The static case 



The Boltzmann equation (2.1) is of course simple to solve in the static case. 
It is easy to verify that 

Po + rA-,{^) (4.1) 



is a constant of the motion. For reasons which will become clear, we now allow 
for a nonzero mass m. Then, in the notation of (3.1), an appropriate solution 



10 



of the Boltzmann equation is 
„ 2Ch 



{2ny 



Sig'^P.Pu - m'/h'))eip,)N [I (po + /MS(x))] , (4.2) 



where N is defined in (2.13). For the Bose case, in order to ensure that F is 
positive, we must restrict A to satisfy 

\rA^{^)\<m/n. (4.3) 



With (4.2), equations (2.5) and (2.6) may be used to define j'^ and T^^. Using 
the change of variable (4.5) below, it follows that j* = 0. With this information, 
one may verify that that equations (1.1) and (1.2) may be integrated to give 



W = —- I d^a; I d^j9 / dJ e{g^"'{x)p^p^ - m^h^ 



{2^y 

xe{po)N[l{po + PAl{^))]. (4.4) 

The p-integration may be done by the substitution 

p'i=Pi- 9ioPo/9oo: (4.5) 

which gives 

9"'p,P. = {9ooryo-9''p[p'v (4-6) 



and thus 



oo 



W = ^J d'^y^ JdlJ dpoH - r'f'N [I (po + r A«(x))] , (4.7) 



where r = rriy/goo/h. 



After expanding A^ as a power series in exponentials, the po-integration may 
be done, to obtain 



oo 

X J2 L^im,/^o/T)n-'[±expi-hrAyT)r, (4. 

n=l 



11 



where 



1 2 

Ln = - {nm^jT^/T) K2{nm^/g^/T) (4.9) 



K2 being the modified Bessel function [8] . These coefficients have the proper- 
ties 



T 

L„ — » 1 for n <^ 



T 

Ln ~ e'X]i{—nra^/g^lT) for n » ^^. (4-10) 



We now consider the Fermi and Bose cases separately. For the Fermi case, 
there is no difficulty in going to the limit m = 0. Then the gluon dependence 
of (4.8) may be expressed in terms of the generalized Riemann zeta function 
$ [81 



dJ ^ n-^ [- e^^{-nPAyT)\' = 

n=l 

fdIexp{-nrA^jT) (!>[-exp{-hPA^jT),4,l]. (4.11) 



The integration over / in (4.11) depends upon the colour group [2]. For the 
case of SU(2), this integration may be made explicit. Writing AqI"" = AqIx, 
we may use 

1 
dI--- = Kb,f f dx---, (4.12) 

-1 



with As = 3/2 and / = \/2 for gluons and A^r = 1 and / = v^/2 for quarks. 
In this case there are only even powers of the gluon field in (4.11). The A^ 
term is the familiar T^ electric gluon mass term. 

One may ask under what conditions the result (4.11) is expected to be relevant 
(aside from the static condition). Since no derivatives of the gluon field appear, 
one expects to need 

k^\rAll k<^T/h, (4.13) 



12 



where A; is a typical gluon wave-number. Also, in order that higher powers of 
Aq than the second should be relevant, we need 

\rA^\ ~ T/h. (4.14) 



Now turn to the Bose case, with condition (4.3). Expanding the exponentials 
in (4.8), we obtain a series of the form 

W = Ch fd^x^ /djf;c„(T,mv/^)(/MS)", (4.15) 



6*00 •' n=0 



where, for m <^ T, 



7r2 /T\^ 1 /T^2 

90 UJ ' ""' "^ 12 U, 



(4.16) 
c?.^ aA-\\og{m./g^lT), c„ -^ a„f-j (m^/^/^)^~" (^ > 3), 



where 03, 04, ■ ■ ■ are numerical constants. For instance, we find that 03 ^ 0.2 
and 04 ~ 0.1. 

The first term in (4.15) corresponds to the action for hot mater in curved 
space-time. The linear term in the gauge field yields a vanishing contribution 
when the integration over colour is performed. The second order term in Aq 
is a good approximation if 

|/M2(x)| < m/n < T/n. (4.17) 



Higher order contributions involving the gluon fields become significant when 
the terms in (4.17) are of comparable orders of magnitude (and (4.13) holds). 
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